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Nuclear Magnetic Resonance (NMR) was successfully employed to test several protocols 
and ideas in Quantum Information Science. In most of these implementations the 
existence of entanglement was ruled out. This fact introduced concerns and questions 
about the quantum nature of such bench tests. In this article we address some issues 
related to the non-classical aspects of NMR systems. We discuss some experiments where 
the quantum aspects of this system are supported by quantum correlations of separable 
states. Such quantumness, beyond the entanglement-separability paradigm, is revealed 
via a departure between the quantum and the classical versions of information theory. 
In this scenario, the concept of quantum discord seems to play an important role. We 
also present an experimental implementation of an analogous of the single-photon Mach- 
Zehnder interferometer employing two nuclear spins to encode the interferometric paths. 
This experiment illustrate how non-classical correlations of separable states may be used 
to simulate quantum dynamics. The results obtained are completely equivalent to the 
optical scenario, where entanglement (between two field modes) may be present. 



1. Introduction 



It is widely accepted that Nuclear Magnetic Resonance (NMR) is an experimental 
technique able to implement an ensemble quantum information processor using 
nuclear spins. From the creation of NMR pseudo-pure states, several quantum 
protocols for quantum information processing (QIP) were implemented in such 
physical system (Oliveira et al, 2007; Jones, 2011). NMR also has turned 
possible the experimental evaluation of abstract theoretical proposals in quantum 
information science as well as tests of principles in foundations of quantum 
mechanics (Auccaise et al, 2012), with a precision rarely obtained through 
other experimental techniques. Thus, NMR has became a powerful experimental 
approach to problems related to QIP, even though it has been criticized due to 



Phil. Trans. R. Soc. A l-[T7l doi: 10.1098/rspa.OOOOOOOO 
This journal is © 2011 The Royal Society 



2 



the absence of entanglement in most of the systems used up to now (OHveira et 
al, 2007), besides the scalabihty issues (Jones, 2011). 

The non-existence of entangled states became a problem when it was 
conjectured the necessity of such a resource to turn possible obtain an advantage 
in QIP (Braunstein et al., 1999). In this sense, as the states used to implement 
quantum protocols in NMR, named the pseudo-pure states, do not present 
entanglement except in singular situations, a questioning about the quantum 
nature of these implementations was raised (Laflamme et al., 2002). It was also 
shown, for a large class of protocols, that the existence of entangled states is a 
necessary (but not a sufficient) condition for quantum computers to present an 
exponential speed-up over the classical ones (Linden and Popescu, 2001), although 
too much entanglement could also be harmful (Gross et al., 2009; Morimae, 2010; 
Tilma et al., 2010). However, there are other important characteristics like the 
effectiveness of the implementation and the manipulation of quantum states, 
which also play some role in the quantum advantage game. Thus, considering that 
the principal aspect of NMR is the excellent control of unitary transformations 
promoted by radiofrequency (RF) pulses, such a spectroscopic technique allows 
one to obtain effective and singular methods in order to manipulate quantum 
states to implement QIP protocols. This is part of the reason of the well succeeded 
NMR-QIP experiments. 

Notwithstanding, as suggested in Refs. (Laflamme et al, 2002; Vedral, 2010; 
Soares-Pinto et al, 2010; Auccaise et al., 2011a; Auccaise et al., 2011b), the 
presence of quantum correlations other than entanglement offers additional 
motivations for the successful NMR-QIP implementations, although, until now, 
not promoting the exponential speed-up of information processing. It is possible 
to measure such non-classical correlations in a bipartite system using a quantity 
called quantum discord (OUivier and Zurek, 2002; Henderson and Vedral, 
2001) [for recent reviews about quantum aspects of correlations beyond the 
entanglement-separability paradigm see Refs. (Celeri et al, 2011; Modi et 
al, 2012)]. One interesting feature of such correlations is that some mixed- 
separable quantum states, or non-cntangicd states, can present non-null discord, 
indicating the existence of non-classical correlations beyond the entanglement- 
separability paradigm. These correlations may have a significant role in quantum 
information science due to the quantum advantage, when compared with its 
classical analogues, obtained by protocols based in non-entangled states (Vedral, 
2010). These concepts allowed us to study experimentally and theoretically the 
existence of such correlations in NMR systems, as well as the effects of the 
decoherence (usually modelled as the amplitude- and phase-damping channels) 
over the system's non-classicahty (Soares-Pinto et al, 2010; Auccaise et al, 2011a; 
Auccaise et al, 2011b; Soares-Pinto et al, 2011). Such a general kind of non- 
classical correlations is a source of quantumness (or a quantum resource) available 
in room-temperature NMR experiments. 

In this article we present an NMR implementation of an analogous of the well- 
known single-photon Mach-Zehnder (MZ) interferometer employing two nuclear 
spins to encode the interferometric paths. MZ interferometer is an important and 
very useful tool for applications in quantum optics as well as for fundamental tests 
of quantum mechanics (Scully and Zubairy, 1997). Such a device enables us to 
determine phase shifts between two paths mediated by two half-silvered mirrors. 
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called beam splitters (BS). At a single-particle level, the wave interference between 
the probability amplitude in the two arms of the interferometer plays a crucial 
role in the phase determination. It is worthwhile to note that it was believed that 
this interferometer, at a single quantum level, needed entanglement (between the 
field modes in both paths of the interferometer) to work properly, in other words, 
for a separable input state of the filed modes after the first beam splitter, the 
output should be also classical. In this case there is no interference between the 
two paths of a single quantum. Experimental implementations of interferometry in 
NMR were already reported in the literature [see for example: (Peng et al, 2003; 
Peng et al, 2005; Auccaise et al, 2012)]. However, such experiments was realized 
using only one nucleus (one spin-1/2 particle) to encode the two interferometric 
paths. Here we present an analogous of a MZ interferometer by means of two 
nuclei that encodes the path information taken by a single quantum. This has 
some advantages, as discussed bellow, but our aim here is mainly to discuss the 
root of the non-classical aspects of NMR systems, as well as to give one more 
indicative of the potential of such a system as a quantum simulator. 

The article is organized as follows. In Sec. [2] we describe how one can quantify 
and measure quantum and classical correlations in highly mixed states in an NMR 
setup. Sec. [3] is dedicated to the detection, but not quantification, of quantum 
correlations. The intention of these two sections is to emphasise the quantum 
nature of NMR systems, since one needs a quantum system to efficiently simulate 
another quantum system. The analogous of the single-photon MZ interferometer 
is implemented in Sec. Finally, in Sec. |5l we present our final discussions. 



2. Measuring quantum and classical correlations in NMR 

For a typical liquid-state NMR system at room temperature, the gap between the 
Zeeman energy levels of the nuclei are much smaller than the average thermal 
energy, e = hwL/2^kBT (with uil being the Larmor frequency), which 

implies that the density matrix can be written as the expansion (Abragam, 1978; 
Oliveira et a/., 2007) 

p^^I + eA/,, (2.1) 

where I is the 2" x 2" identity matrix, n is the number of qubits (encoded in 
spin nuclei) and A/9 is the traceless deviation matrix. This is the well-known 
high-temperature approximation. 

Any manipulation in the above state, like state preparation, quantum state 
tomography, qubit rotations and so on, is performed only over the deviation matrix 
Ap. This is due to the fact that such manipulations are, generally, sequences of 
RF pulses, which are basically unitary transformations [/, that act on the density 
matrix in the following way 

UpU^ = ^l + eUApU^. (2.2) 

Adjusting each pulse length, phase and amplitude, it is possible to obtain a 
very fine control over the density matrix populations and coherences (diagonal 
and off-diagonal elements of the density matrix, respectively, usually defined in 



4 



the computational basis). Together with proper temporal or spatial averaging 
procedures and evolution under spin interactions (Vandersypen and Chuang, 
2004), the RF pulse can be specially designed to prepare any two-qubit 
computational base states, as well as its superpositions, starting from the thermal 
equilibrium state (Bonk et al, 2004; Bonk et al, 2005; Fahmy et al, 2008; Jones, 
2011). Although being possible, in principle, to create any quantum state by means 
of these techniques, it was shown that the obtained state, in general, is a highly 
mixed one that does not possess entanglement since e~ 10~^, thus violating the 
boundary e > 1/(1 + 2^"~^) for the existence of such correlations (Braunstein et 
al, 1999; Linden and Popescu, 2001). 

Nevertheless, such separable state can present other kind of non-classical 
correlations that are conjectured to play a relevant role in QIP protocols (Datta 
et al, 2008; Lanyon et al, 2008; Vedral, 2010). The total correlation contained 
in a state p^s is quantified by the quantum mutual information (Nielsen et al, 
2000; Benenti et al, 2007; Vedral, 2006) 

I{A:B) = S{A) + S{B)-S{AB), (2.3) 

with S'(X) = S'(px) = — Trpx Inpx being the von Neumann entropy and pA = 
Ti bPAB being the reduced density operator for subsystem A, with an 
equivalent definition for partition B. This is a direct generalization of the 
classical mutual information, introduced by Shannon to quantify correlations in 
classical information theory: I{A: B) = H{A) + H{B) - H{AB), with H{A) = 
~'^kPt^^Pk being the Shannon entropy of the probability distribution {p^} 
for the random variable A (Shannon, 1948). In the classical domain, we can 
rewrite this expression in the equivalent form J (A : B) = H(A) — H{A\B), with 
H{A\B) being the knowledge we can get from a random variable A when we have 
measured B, i.e., the conditional entropy. This is the point where the quantum 
and the classical domains breakup. While a classical measure does not disturb 
the system, a quantum one generally does. Therefore, the quantum extension of 
J{A : B) is not straightforward, but a possible one was considered in Ref. (Olliver 
and Zurek, 2002) as 

JiA : B) = SipA) - {pa\b) , (2.4) 

with S^jjB^ {pa\b) = X]j Ij^iP^A) t)6iiig ^ quantum extension of the classical 

conditional entropy 'H{A\B). {H^} is a complete set of projective measurements 

on partition B and p\ = TiB{pAB'^A®^f)/Pi (Pj = Tr[lA CS) H^PAb]) is the 
measured reduced density matrix of partition A. The difference 

V{A : B) = I(A : B) - max J(A:B) (2.5) 
{nf} 

was called quantum discord and is a measure of the quantumness of correlations 
(Olliver and Zurek, 2002, Celeri et al, 2011; Modi et al, 2012), since the 
expressions for X{A : B) and J{A : B) are classically equivalent. We note that 
the second term on the right hand side of Eq. (|2.5p can be regarded as a measure 
of the classical correlations contained in the state pab (Henderson and Vedral, 
2001). 
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Let us back to the NMR system, where all correlations in the quantum state 
PAB comes from the deviation matrix ApAB- In this context, it is desirable to 
express the correlation quantifiers as functions of Ap^B- In order to do this, we 
expand the von Neumann's entropy in powers of the parameter e as 

^(/>) = 2(l-j^TrAp2^+..., (2.6) 

where we have used TrAyO = 0. As the reduced-density operators reads Pa{B) = 
Trij(A)/' = lyl(B)/2 + eAp^(^), with Ap^^^^B) =T^BiA)^PA(B) being the reduced 
deviation matrix, the marginal entropies became 

5(PA(B)) = l-j^TrApi(5) + ... . (2.7) 

Substituting Eqs. (|2.6p and (|2.7p into Eq. (|2.3p and keeping terms up to second 
order we obtain 

^(P)^^ (2Tr (Ap)2 - Tr {ApA? - Tr (Ap^)') , (2.8) 

which is the desired expression for total correlations. 

To quantify the classical correlations we must obtain the measured density 
operator, which is given by (in a symmetric form) (Celeri et al, 2011; Modi et 
a/., 2012) 

= J] (nf nf ) p (nf nf ) 

= ^ + eA7?, (2.9) 

where we defined the measured deviation matrix as Ar] = 
X^jj (j^t {^t ® ^f^' ^^^^ ^bat we applied the projection operators 

on both partitions A and B. This is a symmetric version of the classical 
correlations defined in Eq. (|2.5p and from here on we will adopt such definition 
to quantify correlations [see (Celeri et al., 2011; Modi et al., 2012) for details and 
discussions about non-classical correlation quantifiers]. 

Following the reasoning that led us to Eq. (|2.8p . we obtain the following 
expression for the mutual information of the measured state 

^i^)-^ {2Tr {Al^f - Tr {At^a? - Tr {At^b?] , (2.10) 

and thus for the classical correlation 

C(p)~ max J,(r?), (2.11) 
{nf®nf} 

where At^^^^) = Tr^j-^^) Ar/. The quantum correlation in the composed two-qubit 
system can be directly computed from Eqs. (|2.8p and (|2.1ip as (Soares-Pinto et 
al, 2010) 

Q{p)=X{p)-C{p). (2.12) 
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Figure 1. (Online version in colour.) Relaxation dynamics of correlations. The blue circles are the 
mutual information while the red squares and green triangles represents classical and quantum 
correlations, respectively. In panel (a) we consider two logical qubits encoded in SDS molecule 
under the action of the same phase-damping environment while panel (b) shows the dynamics of 
the two physical qubits of the chloroform. The insets show a schematic diagram of the molecular 
structure of both systems. In both cases we have prepared the same initial state. For details see 
Refs. (Soares-Pinto et al, 2010; Auccaise et al., 2011a). 



Equations (|2.1ip and (|2.12p are general and can be used to quantify quantum 
and classical correlations in every system whose density matrix can be cast in 
the form given by Eq. (|2.ip . Specifically for NMR systems, we have employed 
such expansions to study the relaxation dynamics of correlations in two different 
system: {i) a sodium dodecyl sulfate (SDS) liquid-crystal sample containing ^^Na 
quadrupolar nucleus with spin 1 = 3/2 (Soares-Pinto et al, 2010) and (zi) a liquid 
sample of ^^C enriched chloroform, where we have two spin-1/2 systems (^H and 
^^C nuclei) (Auccaise et al., 2011a). We can use system (i) to encode two logical 
qubits, while system (m) encodes two physical qubits. It is worthwhile to note 
that, during the decoherence dynamics, the amplitude-damping channel (which 
describes energy loss) acts individually on each qubit in both systems, while the 
phase-damping channel (responsible for coherence loss) is a global channel in case 
(i) and it acts independently over each qubit in case (m) (Soares-Pinto et al., 2011; 
Souza et al, 2010). This difference allowed us to detect, in system (ii), a peculiar 
behaviour of the decoherence dynamics of quantum discord, i.e., a sudden change 
in its decay rate (Auccaise et al, 2011a) (theoretically predicted in (Maziero et 
al., 2009)). 

These results reveal two interesting aspects of NMR systems. First, there are 
quantum correlations in these systems, that can be quantified, for instance, by 
the quantum discord. Second, since the relaxation mechanisms present distinct 
features in both systems, they provide a very interesting platform to investigate 
the differences introduced in the dynamics by a global and a local environment. 
An example of this difference is illustrated in Fig. 1 where the relaxation dynamics 
of the correlations between two qubits is shown. While, in the case of a global 
reservoir (Fig. la), we have the expected exponential decay of all kinds of 
correlations (Soares-Pinto et al, 2010), the sudden-change behaviour clearly takes 
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place in the case of independent reservoirs (Auccaise et al, 2011a), as can be 
seen in Fig. lb (see the caption of the figure for details). Moreover, it is quite 
remarkable that the sudden-change still takes place even in the presence of the 
amplitude-damping channel, indicating that such phenomenon could be a strong 
characteristic of the correlations. 



3. Witnessing the quantumness of correlations in NMR 

In the last section we discussed the quantification of quantum and classical 
correlations in a composite systems. This is a very expensive information since 
a full quantum state tomography, beyond the numerical optimization process to 
compute Eq. (|2.1ip . must be performed. However, there are many situations in 
which we do not need to know how much correlations a certain state possess, 
instead it is enough to know only its nature. In order words, we just have to 
distinguish between classical and quantum correlations. In the same spirit of what 
happens to entanglement, some non-classical correlation criteria and observable 
witnesses were proposed in Refs. (Rahimi and Saitoh, 2010; Girolami and Adesso, 
2011; Brodutch and Modi, 2011; Zhang et ai, 2011; Yu et a/., 2011; Maziero and 
Serra, 2010) and experimentally verified in Refs. (Auccaise et al., 2011b; Passante 
et al, 2011; Aguilar et al., 2012). A classicality (or non-classicality) witness is 
regarded as an observable (or a set of observables) that can be directly measured 
in an experimental setup. Depending on its expected value, we know if the state 
has quantum or only classical correlations. The measurement of such a witness 
in an NMR system revealed directly the non-classical nature of its highly mixed 
state (Auccaise et al, 2011b; Passante et al, 2011). 

In what follows, we describe a non-classicality witness and its experimental 
measurement performed in an NMR setup at room temperature (Auccaise et al, 
2011b). In Ref. (Maziero and Serra, 2010) it was shown that the mean value of 
the operator 

3 4 

^p=E E m)p{o,)p\=o, (3.1) 

i=l j=i+l 

is a sufficient condition for classicality of correlations for a wide class of two-qubit 
systems, expressed by 

,0 = — + - Y^iAiaf + af + dat ® af). (3.2) 

i=l 

Here Oi = ® of for i = 1, 2, 3 and O4 = Ei=i(-2if^«^ ^ Wi\^ ® af). af^^^ 
is the zth Pauli operator acting on subsystem A{B). Ai,Bi,Ci, Zi,Wi G JR, with Zi 
and Wi randomly chosen and constrained such that zf = wf = 1. 
In the case of Bell-diagonal class of states 



■11 i \ ^ A R 

Pbd = — + 42^ CiO-i ® (Ti , 



(3.3) 
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^Pbd = is also a necessary condition for the absence of quantumness in the 
correlations of the composite system (for this case (04)pi,^ = (Maziero and Serra, 
2010)). 

One interesting fact about such witness is that it is possible to rewrite the 
observables in Eq. (|3.ip in terms of one component of the magnetization in one 

subsystem as {Oi)p = {a^ I^)?^, with = UA^B[Rni{Oi)pRU{Oi)]U]^^^, where 

Rrii (Gi) = Rn i^i) Rn (^«)) being Rui^^ (Oi) a local rotation by an angle 9i around 
direction rii on subsystem A{B), where 6i = 0, ^2 = ^3 = 7r/2, n2 = y- And = z, 
Ua^b is the CNOT gate with subsystem A being the control qubit. This fact leads 
to a straightforward implementation of this witness in the NMR scenario, since the 
one-qubit magnetizations are the natural observables for these systems. In fact, the 
witness in Eq. (|3.ip was experimentally implemented in a room temperature NMR 
two-qubit system (Auccaise et al, 2011b), directly revealing the non-classical 
aspects of highly mixed states. 



4. Mach-Zehnder interferometer 

An interesting way to test non-classicality in NMR systems can be provided by 
an analogous of the well-known single-photon Mach-Zehnder interferometer. In 
order to perform such an interferometer and test the role of correlations, we 
employ two nuclei to encode the two interferometric path information. As already 
mentioned, this approach differs significantly from the previous interferometric 
measures implemented in the NMR scenario (Auccaise et al, 2012; Peng et al, 
2005; Peng et al, 2003), in which the two-path information are encoded in just 
one nuclear spin. Moreover, due to the fact that distinct nuclei have, generally, 
distinct relaxation times, it turns possible the study of the environment-induced 
phase shift between both paths. This environment-induced phase shift may found 
applications, for example, in thermometry (Stace, 2010) and in the quantum 
illumination protocol (Lloyd, 2008). 

Now, let us briefly review the single-photon MZ interferometer in the optical 
scenario, schematically shown in Fig. 2. In this figure, BS is a 50 : 50 beam splitter 
(one half-silvered mirror), M is a mirror and Da and Db are one- photon sensitive 
detectors. The phase difference (j) between paths A and B can be due to an 
environment-induced phase or just a controlled one. In what follows, the indexes 
A and B must be understood as two spatial field modes (referring to distinct 
paths taken by the photon). 

Let us suppose that the input (pure) state (before the first BS in Fig. 2) is 
given by 

\^,) = \0)a(^\1)b. (4.1) 
The state after the first beam splitter is 

1^-) = i\0)A |1)b + i\l)A <S) \0)b) . (4.2) 

Note that this can be seen as an entangled state of the two field modes (that 
encode path information). The coherent superposition of both paths after the 
first beam splitter is essential for the interferometer to produce an interference 
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Figure 2. (Online version in colour.) Schematic diagram of the Mach-Zehnder interferometer. 



pattern (Scully and Zubairy, 1997). In our experiment, the information about 
both paths will be encoded in two distinct nuclei and, i.e., when talking about 
correlations between these nuclei, we mean the coherent superposition paths A 
and B. 

The final state, after the second BS (see Fig. 2), will be given by 



1^' 



f) 



■ cos 



I ) |0>A^ 



l)B-sin ^ \1)a^\0)b 



(4.3) 



From here on we will omit the tensor product symbol to simplify notation. 
Equation (|4.3p shows the interference pattern between the two paths and reflects 
the wave-like behaviour of the one-photon state. The probability to detect the 
photon in the detector Db (the probability of detect the initial state at the end 
of the interferometer) is given by cos^ {4'/'^)- If ^6 introduce a non-destructive 
detector in any of the paths before the second beam splitter, this interference 
pattern disappears and we observe the particle character of the quantum system. 
This is one the most known representations of Bohr's complementarity principle. 
It is interesting to note that we could modify the present scheme to quantitatively 
study this principle by means of an inequality relating the which- way information 
(particle) and the fringe visibility of the interferometer (wave), as suggested in 
(Englert, 1996). 

We implement such interferometer employing a two-qubit system comprised 
by nuclear spins of and ^^C atoms in a carbon-13 enriched chloroform molecule 
(CHCI3). The sample was prepared by mixing 50 mg of 99% ^'^C-labelled CHCI3 in 
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0.7 ml of 99.8% CDCI3 in a 5 mm NMR tube. Both samples were provided by the 
Cambridge Isotope Laboratories - Inc. The experiments were performed at 25° C 
using a Varian 500 MHz Premium Shielded (^H frequency), at the Brazilian Centre 
for Physics Research. A Varian 5 mm double resonance probe-head equipped with 
a magnetic field gradient coil was used. 

In our NMR version of the MZ interferometer, we have exactly the same 
situation of the optical standard version. The two paths are encoded in two 
different nuclei, with just one quantum of excitation. The nuclear spin in the 
ground state represents the vacuum field state, while the nuclear spin in the 
excited state represents the field excitation (the one photon state) . We choose the 
-^^C nucleus to encode the path in which we apply the controlled phase (path A 
in Fig. 2), while the nucleus encodes the reference path (path B in Fig. 2). 
This choice was motivated by the fact that the Carbon nucleus has the smaller 
transversal relaxation time [see (Auccaise et al, 2011a; Auccaise et al, 2011b) 
for details about the relaxation times for our system]. It is not important for the 
present experiment, but can be explored by experiments studying the effects of 
decoherence in just one path, like the thermometry mentioned above (Stace, 2010). 
Here, the phase difference will be provided by a controlled z-rotation applied to 
the Carbon nucleus. 

The pulse sequence used to implement the analogous of the MZ interferometer 
is shown in Fig. 3. The first block of this figure is used to prepare the deviation 
matrix in a state analogous to |^') = |0)|1), which means that the Hydrogen 
nucleus is in the ground state while Carbon is in its excited state. Some words 
about the meaning of this analogy are in order. In Sec. [2] we have mentioned that 
the spin state is represented by the highly mixed density operator shown in Eq. 
(|2.ip . This happens because the magnetic energy and the thermal energy ratio 
is e~10~^, so the density operator was decomposed in two contributions: the 
white noise term and the deviation matrix term. Recalling that the control over 
the system is given by RF pulses, which are represented by unitary operations 
(besides gradient pulses), we only act in the deviation matrix term (second term 
in Eq. (|2.ip ). We then encode all the information of the state |^') in this matrix. 
At the end we will have a quantum dynamics equivalent to that of a pure state. 

After the state preparation, there are three consecutive blocks that are used 
to implement the MZ interferometer. The first and the third ones play the 
role of the beam splitters. We have constructed a sequence of 128 pulses with 
modulated frequencies and amplitudes by means of the GRAPE optimization 
method (Machnes et al, 2008) to implement the following action on both qubits 



BS = ^ 
V2 



10 

I i 

i 1 

1 



(4.4) 



which represents the beam splitter operation. The time duration of the entire 
sequence in order to implement such operation is 4.5ms. Finally, the three pulses 
between the beam splitters implement a rotation on the Carbon nucleus about 
the z axis by an angle (j). After the pulse sequence depicted in Fig. 3 we performed 
a full quantum state tomography. We repeated the protocol in Fig. 3 for different 
phase shifts. 
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Figure 3. (Online version in colour.) Pulse sequence for the implementation of the analogous of 
a single-photon Mach-Zehnder interferometer in a two nuclei NMR system. The thicker filled 
bars represent 7r/2 pulses, the thinner bars indicate 7r/4 pulses, and the grey bars indicate 7r/6 
pulses with the phases as shown (negative pulse phases are described by a bar over the symbol). 
A = 1/ J and the time periods in the state preparation represent free evolutions under the J 
coupling (Oliveira et al, 2007). 



The experimental results are presented in Fig. 4, where we plot the probability 
of detecting the initial state Pi = \ {^pi\■^p f)\'^ as a function of the phase shift (p. The 
theoretical prediction from Eqs. (|4.ip and (|4.3p are in a good agreement with the 
experimental results showing that the interference expected for the interferometer 
can be obtained in the context of the highly mixed states NMR system. It is 
remarkable that the tiny coherence present in our experiment (of order e) enables 
the interferometer to work properly. The visibility of the interferometer (which 
quantifies the contrast of the interference in any system which has wave-like 
character) is near to unity, without entanglement between the two paths (the two 
nuclear spins in this case). We may argue that the reason for this is the quantum 
correlation of separable states (quantified by the quantum discord) present in the 
composite nuclear spin system. 

The small deviation between the observed and the expected probability in Fig. 
4 is mainly due to the fact that we are continuously applying a RF pulse on the 
Carbon spin to obtain the controlled phase shift between both paths. This pulse 
sequence presents some fluctuations introducing noise into the system. Moreover, 
the pulse is not exactly in resonance with the qubit transition frequency, which 
contributes to the phase mismatch between theory and experiment. We observe 
that there are two decoherence processes involved here. The first one is that caused 
by the usual amplitude- and phase-damping channels, which are always present, 
but that have a negligible effect for the current experiment due to the fact that 
it occurs in a time much shorter than the characteristic relaxation times of the 
system. All the observed decoherence is due to the imperfections of the RF pulses 
employed to generate the phase difference (p between both paths. 

To give some support to our statement that quantum correlations of separable 
states are the responsible for the successful simulation of the single-photon 
interferometer, in Fig. 5 we plot the visibility and the quantum discord between 
both qubits (paths) in a slight modified experiment, where we vary the amount of 
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Figure 4. (Online version in colour.) Probability for the detection of the initial state in the 
experimental implementation of the Mach-Zehnder interferometer in a two nuclei NMR system. 
The (blue) dots are the experimental results and the (red) solid line is the theoretical prediction 
of Eqs. dllj and (gH). 



quantum discord after the first BS operation in the pulse sequence of Fig. 3. To 
obtain this result we performed two experiments. In the first one, we let the system 
to freely evolve (under the natural decoherence) during a certain period of time 
r just after applying the first BS operation and, then, we perform a full quantum 
state tomography. From the obtained deviation matrix we compute the quantum 
discord through Eq. (|2.12p . By varying the time r of the decoherent evolution we 
obtain different amounts of quantum discord, the results are shown by the red 
squares in Fig. 5. In the second experiment we again, after the first beam splitter, 
let the system to freely evolve under the action of the environment. However, 
instead of performing the state tomography, we proceed with the sequence shown 
in Fig. 3 and compute the fringe visibility of the interferometer 

^ ^ max(^/ 1 Z)a I ^/) - min(^/|Z)B|^/) 
max{^f\DA\^f) +mm{^f\DB\'^f)' 

with Da and Db being the detectors operators. As we can see from Fig. 5, the 
decay of the discord is clearly accompanied by the decay in the visibility, showing 
that, without discord, we cannot obtain a visible interference pattern from the 
simulated MZ interferometer. 
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Figure 5. (Online version in colour.) Quantum discord and visibility of the MZ interferometer 
for different decoherence periods r. Note that, although both quantities are represented in the 
same vertical axis, they are completely different quantities. While the visibility is a dimensionless 
quantity, the quantum discord is measures in unities of e^/ In 2 bit. The purpose of the graph is 
only to show that, the fade of the interference pattern is accompanied by the vanishing of the 
quantum discord between both qubits (paths of the interferometer). 



5. Final Discussions 

The QIP either for computing or communication, probably will lead to a 
technological revolution in this century. Having appeared in early 1980s, Quantum 
Computation and Quantum Information evolved faster as an area of basic research 
in physics and as a promise technological alternative. In fundamental physics, the 
concepts of quantum information has stimulated a huge number of new ideas and 
results, which have deepened our knowledge about the nature of the quantum 
behaviour. Particularly exciting are the relations of quantum correlations and 
natural resources for computation and communication. The discovery of quantum 
algorithms, new forms of communication, and cryptography led to a remarkable 
experimental development in order to test these applications. 

Up to now, almost all quantum algorithms and communication protocols have 
been tested on systems with a few qubits. However, for large-scale applications 
of quantum devices it is necessary to study the generation, handling and the 
storage of quantum states in such systems. We may say that the experimental 
technique that have stood out in this context was NMR, because it allows a quite 
precise manipulation of quantum states. This precise control of spins states is due 
to the radio-frequency technology developed over decades. It has been used in 
an inventive way for NMR methods and contributed to the impressive and fast 
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success of the technique in QIP (Ohveira et al, 2007; Jones, 2011). NMR also 
allows, through the manipulation of nuclear spin qubits, the simulation of more 
complex systems (Alvarez and Suter, 2010; Du et al, 2010; Lu et al, 2011; Souza 
et al., 2008; Souza et al, 2011). 

Although all the successful use of NMR in QIP, some years ago it was shown 
the impossibility of existence of entangled states in such experimental technique 
at room temperature. By that time, it was supposed that this kind of quantum 
correlation was intrinsically linked with the speed-up of the quantum algorithms. 
Thus, one question could be raised: What is the quantum aspect of NMR states 
that allows quantum dynamics and its application in QIP? This question is 
partially answered by the existence of quantum correlations in separable states, 
as revealed, for example, by the quantum discord. As we have argued in this 
article, such kind of correlations seems to occur naturally in room temperature 
NMR systems, as far as it can be easily quantified or witnessed by the methods 
discussed here, and it is responsible for the quantumness of the system. The 
analogous single-photon MZ interferometer presented in the last section provides 
a clear illustration of the non-classicality present in NMR highly mixed states. The 
coherence present in the deviation matrix allows for a path interference analogous 
to that one obtained in the optical context where an entangled state (between 
two field modes) is present. In other words, we may argue that the separable 
state (non-entangled) quantum correlations present in the two qubit NMR state 
are the resource that enables the interferometer to work properly with a visibility 
near to unity. 

In summary, the phase coherence between two-nuclear spins in the NMR 
highly mixed states can encompass non-classical correlations (as measured by 
the quantum discord). Such quantum correlation are available for performing 
quantum dynamics (path interference) allowing the simulation of several systems, 
as for instance, the single-photon MZ interferometer. The discussions presented 
here may opens the way for a new and very exciting tests of the quantum aspects 
of nature and of protocols in information science that exploit this kind of non- 
classicality, which is present even in the NMR highly mixed states scenario. 
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